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ABSTRACT 

We present millimeter interferometric observations of the molecular envelope around the red super- 
giant VY CMa with the SubMiUimeter Array (SMA)B The high angular resolution (< 2") allows us 
to derive the structure of the envelope as observed in the 1.3 mm continuum, ^^C0(2-l), ^'^00(2-1) 
and 80(65-54) lines emission. The circumstellar envelope is resolved into three components: a dense, 
compact and dusty central component, embedded in a more diffuse and extended envelope plus a 
high velocity component. We construct a simple model, consisting of a spherically symmetric slowly 
expanding envelope and bipolar outflows with a wide opening angle 120°) viewed close to the line 
of sight (i — 15°). Our model can explain the main features of the SMA data and previous single-dish 
CO multi-line observations. An episode of enhanced mass loss along the bipolar direction is inferred 
from our modelling. The SMA data provide a better understanding of the complicated morphology 
seen in the optical/IR high resolution observations. 

Subject headings: circumstellar matter - supergiants - stars: individual (VY Canis Majoris) - radio 
lines: stars 



1. INTRODUCTION 

The terminal stage in the evolution of massive stars (M 

> 10 M0) is characterized by very short timescales (~ 
10** yr) and drastic changes in their immediate circum- 
stellar environment. In the red supergiant phase, periods 
of large mass loss release a thick circumstellar envelope 
of molecular gas and dust. Understanding the mass loss 
process from red supergiants is extremely important in 
that it will eventually govern their evolution as well as 
the structure of their surrounding envelope. In a galac- 
tic context, supergiants contribute to the enrichment of 
the interstellar medium in dust and heavy elements and 
participate to the global galactic chemical evolution. 

VY Canis Majoris (VY CMa) is one of the brightest 
known red supergiants. Its distance has long been de- 
bated (Hcrbig 1970; Lada & Rcid 1978), but there is now 
a consensus on a value of 1.5 kpc, consist ent with the ve- 
locity and proper motions of H2O masers ([Richards et al.l 
I1998I). The stellar lum inosity is estimated to ~ 2 10^ L© 
( Monnier et al.|[l999r ). The stellar mass is thought to be 

> 15 M0. The basic properties of VY CMa are summa- 
rized in Table 1. 

VY CMa is characterized by an apparent very high 
mas s loss rate, of order 2-4 10~^ M^/yr (see 
e.g. 'Bowers et al."1983': [Panchi et al.lll994l : iStanek et all 
[1995; Harwit ct al. 20Qil), which produced a very thick 
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circumstellar envelope of dust and molecules, which 
almost completely obscurs the central star. The 
optical/infrared morpholgy of the nebulosity around 
VY CMa is very complex, showing arcs, fila ments, bright 
knots in a very complicate d arrangement (jSmith et al.l 
I2OOII: iMonnier et al.l I1999D. Long slit spectroscopy 
(|Smithll2004l : lHumphrevs et al.ll2005D reveals a very com- 
plex kinematics w ith evidence for lo c alized and episodic 
mass loss events. [Humphreys et all (|2005[ ) showed that 
the prominent arcs are expanding relative to the star and 
represent discrete ejection events. The ejection dates 
are ranging from 200 to 1000 years ago and state that 
VY CMa has had a period of major activity during the 
past 1000 years. These works however rely on scattered 
lights which may reflect only the illuminated part of the 
nebula and the origin of the features seen in the optical 
and near infrared images is very uncertain. 

We note that in a recent paper. E3ecin et aTl (l2006f ) pro- 
pose a theoretical model, based on a spherically symmet- 
ric envelope, that requires time varying mass-loss rates 
to reproduce previous single-dish CO multi-line spectra. 
The high angular resolution view offered by the SubMil- 
limeter Array (SMA) now provides new constraints on 
the structure and kinematics of the circumstellar molec- 
ular envelope of this peculiar object. In this paper, we 
present our SMA observations and discuss about the 
morphology and kinematics of the molecular envelope 
aroimd VY CMa. We propose a simple model to explain 
our data and compare with observations at other wave- 
lengths. 

2. OBSERVATIONS 
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We obtained 1.3 mm interferometric obser vations of 
VYCMa with the SubMillimeter Array (SMA. lHo et alJ 
|2004|) . which consists of eight 6 m antennas, located on 
top of the Mauna Kea, Hawaii. The phase reference of 
the array was set to R.A. = 07''22™58!328 and Dec. = 
-25°46'03'.'19. 

VY CMa was observed at 2 different epochs and with 
2 different configurations. The first observation was car- 
ried out on 2004 July 31^*", with 7 antennas in a com- 
pact configuration, including projected baselines from 8 
to 99 m. The system temperatures at 1 mm were be- 
tween 100 and 200 K for the different antennas. In order 
to improve the angular resolution, a second track was 
observed on 2005 January 19*^, with 7 antennas in an 
extended configuration. The projected baselines ranged 
from 18 to 192 m. The system temperatures were be- 
tween 200 and 400 K. The bandpass was calibrated using 
Uranus and/or Jupiter, while Callisto set the flux cali- 
bration. The nearby quasars 0607—157 and 0730—116 
were regularly observed for the gain calibration. After 
calibration, the visibilities from both tracks have been 
combined together. The data reduction and calibration 
were done under MIR/IDL. Imaging and data analysis 
were entirely processed under GILD AS ^. 

The heterodyne receivers were tuned to allow multi- 
line observations including the CO (2-1) line in the up- 
per side band (USB) and the 1^00(2-1) and SO(65-54) 
lines in the lower side band (LSB). The SMA correlator 
bandwidth is 2 GHz wide for each sidebands, which are 
separated by 10 GHz. The instrumental spectral resolu- 
tion was 0.8125 MHz, but we smoothed it to 3.25 MHz 
{i.e. ~ 4 km s-i) for the 1^00(2-1) and SO(65-54) lines 
and to 8.125 MHz {i.e. ~ 11 km s~^) in the case of 
i3CO(2-l). 

We obtained the map of the 1.3 mm (225 GHz) con- 
tinuum emission by averaging the line- free channels from 
both LSB and USB, resulting in a total bandwidth of 3 
GHz. To maximize the angular resolution, we applied 
uniform weighting to the continuum visibilities, leading 
to a synthesized beam size of 1.38" x 1.23" (PA = 169°). 
We estimate a rms noise level of 8 mJy/beam in this 
continuum map. 

The channel maps of the 1^00(2-1) and SO(65-54) 
line emission were obtained using natural weighting 
and gaussian tapering with a full width at half max- 
imum (FWHM) of 130 m, making a compromise be- 
tween brigthness sensitivity and angular resolution. The 
^■^00(2-1) map was derived using a taper distance of 80 
m to increase the signal-to- noise ratio. The resulting 
synthesized beam sizes and la noise levels are reported 
in Table 2. The continuum visibilities were subtracted 
before imaging. All the maps were deconvolved with a 
standard CLEAN/CLARK algorithm. 

At 230 GHz, the FWHM of the 6 m SMA antennas 
primary beam is 54". The shortest projected baseline 
being 8 m, we may loose informations on structures larger 
than 40". Nevertheless, the resulting coverage of the 
uv-plane is rather scarce until 18 m. We may therefore 
also miss part of the emission associated with structures 
larger than 18", although we show in i )3.2.2l from the 
comparison of the ^^C0(2-l) SMA and single dish data, 
that the missing flux is limited to less than 20 %. 

1 http://www.iram.fr/IRAMFR/GILDAS/ 



3. RESULTS 
3.1. l.S mm continuum emission 



iLipscy et al.l (|2005l ) observed VY CMa at centimeter 
wavelengths using the VLA. They do not resolve the 
source at an angular resolution of ^ 1" at 1.3 cm. They 
also present Keck mid-infrared imaging and show that 
the dust emission is elongated in the south-west direction 
(PA = 227°) with FWHM ~ 0.4" and 0.3" at 11.7 and 
17.9 /im respe ctively. From diff r action -limited speckle 
interferometry, IWittkowski et all (jl998[ ) measure a dust 
shell size of 138 x 205 mas (2.17 /xm), 80 x 116 mas 
(1.28 /xm) and 67 x 83 mas (0.8 /im) with roughly a 
north-south elongation. Thus, at near-IR wavelengths, 
the size of the dust-emitting region decreases towards 
shorter wavelengths, progressively tracing warmer dust 
close to the star. 

As seen from Fig[Tl the continuum emission is not re- 
solved by our synthesized beam (~ 1.4"). The continuum 
visibilities are well fitted by a point source of 288 ± 25 
mJy (with 3<t rms), located at a position -1-0.15" ± 0.01" 
to the east and —0.05" ± 0.01" to the south, relatively to 
the phase center. Our measurement of the flux density is 
consistent with the value of 27 zt 40 mJy from prev ious 
SMA observations at 215 GHz (I Shinnaga et al.ll2004[ ) but 
differs from the previous bolome tric values of 630 zb 65 
mJy obtained with the JCMT (Knapp et al."1993) and 
378 ± 20 mJy with the SEST (Walmsley 1991). Part 
of the difference could be explained by contamination 
from line emission in the bolometric observations. We 
note, however, that even the ^^C0(2-l) line, presumably 
the strongest in this frequency band, contributes only 
for about 2% and 7% of the total flux density measure- 
ments with the JCMT and SEST, respectively. On the 
other hand, we cannot exclude the presence of an ex- 
tended low surface brigthness emission that could have 
been missed in our observations. Independant checks of 
the continuum flux density in the upper and lower side 
bands give consistent values within 10%, implying that 
our measurements are not severly contaminated by weak 
line emission. 

Because the continuum emission is of order of ^ 1 
mJy in the centimeter wavelenghts (jLipscv et al.ll2005l) 
and given that the contribution of the central star (black 
body emission) amounts for S, ~ 10 mJy, i.e. is negligi- 
ble with respect to the total 1.3 mm continuum emission 
S^, this latter should arise almost completely from dust 
thermal emission. Under the Rayleigh- Jeans approxima- 
tion and if the dust emission is optically thin, we can 
estimate the dust mass M^ within our 1.4" beam, ac- 
cording to: 



2c^D^agPg{S,-S,) 



(1) 



where D is the distance to the source, pg and a^ are the 
mass density and radius of the dust grains, is the 
grain emissivity and the average temperature of the 
dust envelope. We assume typical values for oxygen rich 
stars, dug = 0.2 pm, pg = 3.5 g.cm~'^ , = 5.65 10~*(i^ / 
274.6 GHz), following [Kl^appet al.l (|1993f ). The ratio of 
IR flu xes indicates a tem perature T^ ~ 250 K in the dust 
shell (jLipscv et al.ll2005f) . Accordingly, we derive a dust 
mass of 1.5 1 0"^ M^, close to the value M^ = 1.3 10"^ 
M0 given by iLipscv et al.l (|2005f ). Assuming that the 
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dust is uniformly distributed within a radius comparable 
to our synthesized beam (upper limit of 1.4"), the optical 
depth towards the central star would be ti.3,„„j ^ 0.01. 

Furthermore, assuming an expanding velocity Vcxp — 
15 km s^^ (similar to the escape velocity of ~ 20 km s~^ 
at the dust formation radius, see also S|4|) leads to an age 
of order of 300 years. The average mass loss rate over 
this timescale is 4.7 lO"** M0/yr, assuming an isotropic 
mass loss and a dust-to-gas ratio of 100. This value is 
slightly hi gher than previous measurements for VY CMa 



(see e.g. lOanchi et all Il994l : iHarwit et ahl l2001h and 



comparabl e to mass loss rates es ti mated for ot h er su - 
pergiants_(!Netzer & Knaod llMZl) • ISmith et all (|2001l ) 
and iHurrrph rcvs et al. (2G0ft) showed that localized and 
short-term ejection events in VY CMa could in fact reach 
10-3 Me/yr. 

On the other hand, as pointed out by the referee, there 
is a possibility that the continuum dust emission origi- 
nates from a circum stell ar disk, as advocated earlier by 
iBowers et all (pSS) and Ri chards et all (|1998f ) from ob- 
servations of OH and H2O maser emissions. However, 
the angular resolution of our current data is not enough 
to constrain this possibility. 

3.2. Molecular line emission 
3.2.1. Description of the molecular line emission 

The 1^00(2-1) emission (Figij) extends to - 7" 
around the central star and is detected over a (LSR) 
velocity range —25 to 70 km s~^. There are three dis- 
tinct emission peaks in the channel maps: a first one at 
a velocity of about —5 km s^^ with an offset of 1" to 
the east, a second one close to the systemic velocity of 
VY CMa (VsYS 22 km s^^), centered on the stellar 
position and a third one around -f45 km with an off- 
set of 1" to the west. The line emission is stronger in the 
western redshifted part of the nebula, where it reaches ~ 
5 Jy/beam. The velocity channels around the systemic 
velocity of VY CMa (12 < Ylsr < 33 km s^^) show a 
north-south elongated structure of roughly 5" x 3", em- 
bedded in a more diffuse envelope. A second weak peak 
is seen 3" to the west. 

Our ^^00(2-1) channel maps (Fig|3]), although of low 
signal-to-noise, present globally the same features as seen 
in ^^C0(2-l). The two peaks are well detected at veloci- 
ties ~ —5 and -1-50 km s~^, and separated by ~ 2". The 
emission around the systemic velocity is weakly detected. 
The channel maps of the 80(65-54) emission (Fig|4]) also 
show the same features as that of the ^^C0(2-l). The 
integrated intensity maps for the three species are pre- 
sented in FiglH 

The morphology and kinematics of the ^^C0(2-l) emis- 
sion around VY CMa do not reproduce those of a typical 
spherical expanding shell. We interpret the east and west 
offset bright peaks, with high velocity relatively to the 
systemic ve locity, as the signatur e of a east-west bipo- 
lar outflow. IShinnaga et al.l (|200^ already suggested the 
presence of a bipolar outflow around VY CMa from the 
butterfly like shape of the SiO(l-O) emission. The idea 
was further suppor ted by SiO masers li near polariza- 
tion measurements (jShinnaga et al.ll2004f ). which result 
in a mean polarization position angle of 72° , in the axis 
of t he SiO bipolar structure. Previous observations of 
OH ([Bowers et al.lll983lll993l) and H2O (|Richards et all 



IT998I) masers also show some elongation in a similar di- 
rection with position angles of 50° and 60°, respectively. 
However, the maser emission is clumpy and might not 
trace the complete geometry of the outflow. 

In addition, the 1.3 mm continuum and CO emission 
around the systemic velocity indicate the presence of 
a central dense component. The position- velocity dia- 
grams (Fig|6]) , along the east- west and north-south direc- 
tions, i.e. along and perpendicular to the outflow axis, 
support our interpretation. Along the east-west direc- 
tion, the high velocity bipolar outflow and the slowly ex- 
panding shell appear clearly as two distinct kinematic 
components. The bipolar outflow emission shows a 
point symmetry relative to the stellar position. It ex- 
tends to roughly 5" from the star and has a velocity 
up to 45 km s^^ with respect to the systemic veloc- 
ity. This corresponds to a timescale of about 800 yr, 
compa rable to the past period of strong mass l oss ac- 
tivity (jHumphrevs et al.l 120051 : iDecin et al.ll2006l ). The 
slowly expanding envelope can be seen out to a radius of 
roughly 5". Given an expansion velocity of 15 km 
(see 21), the corresponding timescale is 2400 yr. The 
overall features are similar in the ^^C0(2-l) and 80(65- 
54) position- velocity diagrams, indicating that they arise 
from the same regions. 

3.2.2. Comparison with single dish spectra 

We present in FigjT] the spectra of the 1^00(2-1), 
^•^C0(2-l) and 80(65-54) lines, taken at the central posi- 
tion and convolved to a beam size of FWHM 19'.'7, equal 
to the JCMT beam at 230 GHz. The 1^00(2-1) line 
profile obtained from our interferometric data can then 
be dir ectly compared to th e single-dish -^^C0(2-l) spec- 
trum (|Kemper et al.l 120031 ) at the same angular resolu- 
tion. The CO (2-1) 8MA spectrum shows three peaks 
at Ylsr -5, 20 and 45 km s-\ simflar with the JCMT 
spectrum. The high velocity channels {Vlsr < 5 and > 
35 km s""'^) reach the same brightness temperature as 
the single-dish observations, indicating that the interfer- 
ometer recovered all of the emission for these velocities. 
Close to the systemic velocity of VY CMa, i. e. for 5 < 
ViSfl < 30 km s~^, there is a small difference between 
the two profiles, although less than ^ 20% of the emis- 
sion is lost by the interferometer. We estimate the inte- 
grated intensity of the ^^C0(2-l) line to 65 K km s"^, 
co mparable to th e value of 66 ± 7 K km reported 
by iKemper et al.l (j2003D . The integrated intensities of 
the ^^00(2-1) and 80(65-54) fines are estimated to 4 
K km s~^ and 28 K km s~^, respectively. The 3cr rms 
uncertainty is 1 K km s~^. 

3.2.3. SO emission 

While carbon chemistry is mostly absent in the circum- 
stellar envelope (C8E) of 0-rich stars, most of the carbon 
being locked up in CO, sulfur chemistry is expected to 
be relatively active. 8-bearing species are indeed com- 
monly observed in the massive C8E around 0-rich stars 
(iSahai fc Wanniei][l99llOmont et al.l[T993h . 

Und er chemical thermodynamical equilibrium, iTsujii 
(jl973() showed that the sulfur in C8Es is mostly 
in the f orm of H28. La t er non-equilibrium chemical 
models (8calo & 81avskJ \imi iNeiad fc Millml [19881 : 
iWillacv fc Millar ,1997. ) predict that a large proportion 
of sulfur should go into SO and SO2, as a result of the 
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reactions; S + OH — > SO + H, HS + O — > SO + H and 
SO + OH — > SO2 + H. OH is a natural product of the 
photodissociation of H2O, and VY C Ma exhibits intense 
maser emission for th ese two species ([Bowers et al.|[l983t 
[Richards et alJ l998). On the other hand, the production 
of HS requires highly endothermic reactions. From our 
SMA observations, the SO and CO morphology and kine- 
matics are similar, implying that SO is present both close 
to the star and in the b i polar outflow. Observations of 
ISanchez-Contreras et al.l (|200 J) towards the bipolar neb- 
ula OH 231.8-1-4.2 show a similar trend. These authors 
argue for the presence of SO both in the central com- 
ponent and in the accelerated lobes, with no significant 
difference in abundance between the two components. 

Assuming optically thin emission and an excita tion 
temperature of T^^ = 40 K (Sahai fc Wannieilll992f) . we 
estimate a SO column density ~ 10^^ cm~^ towards the 
center. The corresponding SO abundance relative to H2 
is - 10-^. 

High angular resolution observations of other S-bearing 
species such as H2S and SO2 would provide a better un- 
derstanding of the sulfur chemistry and put stronger con- 
straints for chemical models. 

4. A SIMPLE MODEL OF THE ENVELOPE 

We have constructed a simple model in order to better 
understand the geometry and physical conditions within 
the circumstellar envelope around VY CMa. The com- 
bination of our high angular resolution C0(2-l) data to- 
gether with previous singl e dish observati ons of multiple 
CO rotational transitions (jKemper et al.l i2003) provides 
very useful constraints for our model and leads us to a 
coherent picture of the envelope. 

The geometry of our model is illustrated in FiglHl The 
envelope is assumed to be axisymmetric and radially ex- 
panding. The bipolar outflow direction defines the sym- 
metry axis of the envelope. We project the envelope di- 
rectly into a regular 3-D grid. At each grid point, the 
gas density, kinetic temperature and expansion velocity 
can be specified. The emerging CO intensity along each 
ray through the envelope is calculated by integrating di- 
rectly the standard radiation transfer equation. We take 
into account 15 rotational levels of the CO molecule in its 
vibrational ground state. The populations on the differ- 
ent rotational levels, which are necessary for the calcula- 
tion of the line opacity and source function at each grid 
point, are determined by solving the statistical equilib- 
rium equations within the framework of the large velocity 
gradient formalism. The collision rate s between CO and 
molecular hydrogren are taken from iFlower fc Launavl 
11985) and calcul ated for different temperatures using 
the prescription of lde Jong. Chu fc Dalgarnol ()1975[ ). Be- 
cause the envelope around VY CMa is oxygen rich, we 
adopt an abun dance f = [C0 ] /[H2] = 3 10"'' as used in 
the analysis of iKemper et al.l (|2003f l. The local line pro- 
file, which is needed to integrate the radiative transfer 
equation, is determined through the turbulence velocity 
of the molecular gas. We assume a turbulence velocity of 
1 km s~^ in the slowly expanding envelope and 3 km 
in the bipolar lobes, respectively. This assumption is 
reasonable as a higher outflow velocity in the lobe would 
result in more turbulent environment. 

As discussed in i j3.2.11 two distinct kinematic compo- 
nents are clearly present in the nebula around VY CMa: 



(1) a slowly expanding envelope slightly elongated in the 
north-south direction surrounded by a diffuse halo with 
some extension toward the north-west direction and (2) 
a high velocity bipolar outflow oriented in the east-west 
direction. However, to keep our model as simple as pos- 
sible, we assume that the slowly expanding component is 
spherically symmetric. We adopt an expansion velocity 
of 15 km s~^ for this component, which allows us to sat- 
isfactorily fit the central peak seen in single-dish CO ob- 
servations. The presence of diffuse emission suggests that 
the gas density in the slowly expanding envelope peaks 
more sharply toward the center than for the case of con- 
stant mass loss. We therefore adopt a simple density dis- 
tribution n(r) — 4.5 10^ (10^^ cm/r)^-^, which deviates 
slightly from the usual density distribution generated 
by constant mass loss. This adopted distribution corre- 
sponds to a mass loss rate M — 4.5 10"^ (10^^ cm/r)°-^ 
Mq /yr, or a decrease in mass loss rate by a factor of 4 be- 
tween the inner radius (Rin and the outer radius (Rout)- 
The kinetic temperature of the molecular gas in the 
oxygen-rich circumstellar envelopes with constant mass 
loss has been known to roughly foll ow a power law T(r) 
~ 7-~" fro m the modelling w o rks of iGoldreich fc Scovilld 
(|1976f) and lJusttanont et al.l (|1994f) . Although the mass 
loss rate in the slowly expanding envelope of VY CMa 
is not strictly constant, resulting probably in some slight 
change in the temperature profile in comparison to the 
simple constant mass loss case, we expect that qualita- 
tively, the temperature profile still follows the power law. 
Therefore, we assume in our calculations, for the sake of 
simplicity, that the kinetic temperature varies as T(r) — 
50 K (2.25 lO^'^ c m/r)" '^, which is similar to that used 
by IKemper et all (|2003D to model the oxygen-rich cir- 
cumstellar envelopes. 

In order to reproduce the velocity gradient seen in the 
position- velocity diagrams (Figl6|), the outflow velocity 
in the bipolar lobes is allowed to increase linearly with 
radial distance from the central star. In addition, the 
outflow velocity in the bipolar lobes also varies as a func- 
tion of the polar angle 9 measured from the bipolar axis. 
The strong emission peaks at large velocity shifts sug- 
gest the presence of an enhanced density shell within 
the bipolar lobes. However, this density enhancement is 
more pronounced in the receeding lobe as clearly seen 
in the position-velocity diagram. With the density en- 
hancement included in the model, we find that the bipo- 
lar lobes and the slowly expanding shell can account for 
the highly unusu al triple peak profile of single-dish CO 
lines observed bv IKemper et al.l ()2003f ) (see Fig|9|). 

The intensities of the CO lines calculated along indi- 
vidual lines of sight are used to form the model channel 
maps fFig fTTj) . We convolve these model channel maps 
with a gaussian beam with specified FWHM to produce 
the model CO line profiles, which are indicated by the 
dash ed lines and compared with the JCMT observations 
of K emper et a l. (2003) in FigjS 

We infer from our modelling that the velocity of the 
gas in the lobes increases linearly from 15 km s~^, which 
is similar to the expansion velocity of the spherical en- 
velope, to 45 km at the outer radius. The bipolar 
outfiow is found to be very open, with an opening angle 
of 120°, and inclines at 15° from the line of sight. All 
the model parameters, which provide a reasonable fit to 



The molecular envelope of VY CMa 



5 



both the SMA and single dish data, are shown in Table 
3. 

Although relatively simple, our model provides a qual- 
ita tive understanding o f the single dish data obtained 
by iKemper et al.' (2003") and also of the structure of the 
envelope (FigslH HH The high angular resolution 

SMA data show that the C0(2-l) emission originating 
from the high velocity bipolar lobes produces a double 
horn profile at roughly ±20 km s~^ from the systemic 
velocity. The central broad peak of the C0(2-l) line 
comes from the slow component. Higher transitions of 
CO probe progressively warmer and denser gas located 
in the inner part of the envelope. However, both the den- 
sity and temperature of the gas in the inner region of the 
bipolar outflow, as suggested by our model, are lower 
than that in the slowly expanding envelope. As a re- 
sult, the central peak becomes more prominent in higher 
transitions while the strength of the emission from the 
bipolar lobes decreases noticeably. In addition, the lower 
outflow velocity in the inner part of the bipolar lobes also 
leads to a narrower full width at zero velocity for higher 
transitions of CO. This behavior can be clearly seen in 
Figil 

By comparing the model prediction for the strength 
of the ^'^C0(2-l) line with our SMA observations, we 
estimate the isotopic ratio ^^C/^'^C — 60. This ratio 
is significantly larger than the values ^^C/^'^C = 6 and 
12 found i n the red supergiants a Ori and a Sco, re - 
spectively (jHarris fc LambertJ[l98l iHinkle et al.l[l97l . 
It is likely that, prior to the current mass loss episode, 
VY CMa experienced no significant mixing of H burning 
products to the surface. Our model predicts that most of 
the ^•^00(2-1) line emission should arise from the bipo- 
lar outflow fFig fTU]) . which is consistent with our SMA 
imaging data (Figlj). 

We emphasize here that our model is constructed to 
explain the general features of the high resolution SMA 
data and single-dish CO observations. Other features 
such as the slight elongation of the slowly expanding shell 
and the more complicated structure within the bipolar 
lobes are not considered. However, our model could serve 
as a useful starting point when higher angular resolution 
data of CO (2-1) and higher transitions become available 
in the future. 

5. DISCUSSION 

5.1. Connection with optical imaging /spectroscopy data 

Our high angular resolution C0(2-l) data bring a 
clearer overall picture of the molecular envelope around 
the enigmatic red supergiant VY CMa. The envelope 
consists of two kinematic components: a slowly expand- 
ing shell slightly elongated in the north-south direction 
and a high velocity bipolar outflow oriented in the east- 
west direction with small inclination angle to the line of 
sight. The presence of these two kinematic components 
have already be en hinted in previous observations of H2 O 
and OH masers (|Bowers et al.lll983l[l99llRichards et al.l 
[l998). Given the small inclination angle ~ 15° of the 
bipolar outflows to the line of sight, as derived from our 
model and assuming uniform illumination in the neb- 
ula, one would expect that the envelope appears only 
slightly asymmetric with a brighter east lobe (nearer 
side). However , optic al a nd near IR linages o btained by 
iMonnier et"al1 (|1999f) and lSmith et all ()2001[ ) aU show a 



highly irregular and asymmetric nebula around VY CMa. 
The nebula is seen to the west and southwest of the bright 
point-like source identified as the obscured central star 
VY CMa and is likely the result of a non-isotropic illumi- 
nation by the central star. Comparing the optical images 
with our SMA data, we note that the optical nebulosity 
corresponds to the location of the west (far side) lobe 
and just a small part of the east (near side) lobe. It 
is therefore likely that the light from VY CMa escapes 
preferentially in the direction of the west bipolar lobe 
and scattering on dust brings a fraction of the light back 
toward the observer. 

The strong 1.3 mm continuum emission detected in our 
data (see H3.ip confirms the presence of a thick inner dust 
shell surrounding VY CMa. The average optical depth 
would be so high, if the distribution of material in the 
dust shell is smooth, that all the stellar light would be 
completely extincted. The complex nebula seen in the 
optical and near IR then suggests that stellar light must 
be able to escape through some opening or holes in the 
dust shell, creating search light beams that illuminate the 
nebula in an irregular pa ttern. We note th at a similar 
suggestion was made by IWallersteinI (|1978f l to explain 
the subtle but peculiar change in the relative position 
between the star and the brightest knot in the nebula. 

The thick and compact inner dust shell detected in 1.3 
mm continuum emission might be the region wh ere TiO 
and ScO bands (|Herbiglll974t IWallersteinI [l986l) seen in 
emission originate. The calculated excitation tempera- 
tures of these molecular bands range from 600 to 800 K, 
suggesting that they come from a hot and dense region 
located close to the star. 

High resolu t ion o ptical and near IR im ages of 
IMonnier et all (|1999( ) and ISmith et all (|2001l ) show a 
very prominent long curved tail or arc to the north- 
west of VY CMa together with several fainter arcs and 
condensations in the envelope. The location of the long 
northwest arc roughly coincides with the extended and 
somewhat diffuse C0(2-l) emission in the western part of 
the envelope (Figl2|) . Together with strong central emis- 
sion core, this extended emission is conspicuous in all 
channel maps redshifted with respect to the central star, 
most prominent in the velocity range 21 - 33 km and 
46 - 54 km s^^. The orientation of this feature changes 
from west to north-west with increasing velocity upon 
closer inspection. We can identify a similar feature in 
the channel maps of the SO emission at the same veloc- 
ity, although there are some slight differences at higher 
redshifted velocities, which are probably related to local 
excitation of the SO molecule. It is therefore possible 
that the scattering of stellar light by dust in this struc- 
ture produces the prominent northwest arc seen in opti- 
cal and near IR images. Future observations of molecular 
lines at higher angular resolution and detailed modelling 
of the light scattering by dust inside the envelope will 
help to clarify the relation between the material distri- 
bution and the optical and near IR appearance of the 
nebula. 

Recently, iHumphrevs et al.l ()2005[ ) have obtained high 
resolution long slit spectroscopic data of VY CMa. By 
measuring the velocity of the refiected stellar absorption 
lines and K I emission lines they suggest that the promi- 
nent northwest arc is expanding at a velocity of ~ 50 
km s^^ on the plane of the sky. More surprising is the 
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existence of diffuse gas and dust more or less stationary 
with respect to the central star, mostly along the slit 
positioned to the southeast of VY CMa (slit V in their 
Fig.l). Thus, components with very different kinematics 
seem to co-exist in the nebula. To reconcile the optical 
spectroscopic data with our SMA results, we note that 
for the simple case of single scattering, the light scattered 
by moving reflective material (dust particles) within the 
envelope is Doppler shifted by Voxp(l — cosw), where 
Vcxp is the expansion velocity and lo is the angle mea- 
sured from the line of sight to the radial vector of the 
scattering point (see Fig[5]). The velocity shift can be 
seen to be largest in the west lobe, which is receeding 
along the hue of sight, because w > 90°. If the associa- 
tion of the extended feature in the west lobe as traced by 
CO emission and the northwest arc is correct, we would 
expect the velocity shift induced by scattering of mov- 
ing dust to vary from 15 km s~^ for material in the arc 
moving close to the plane of the sky, which is part of 
the slowly expanding shell, to ~ 50 km s~^ or higher for 
material within the bipolar lobe. This simple estimate 
of the velocity shift i s cons istent with the measurement 
of iHumphrevs et all ()2005l ) . For absorption lines seen 
in scattering through the east lobe, the velocity shift 
is expected to be close to zero because the angle u is 
small. That might explain the presence of dust and gas 
al most stationary wi t h resp ect to the star as suggested 
bv IHumphrevs et al.l (|2005( ). However, similar behavior 
seen in resonant K I lines are still difficult to understand. 

5.2. Mass loss history 

There are strong evidences that the mass loss rate 
in VY CMa is time variable, such as the prom inent 
arc- like features seen in optical (jSmith et al.ll200ll ) and 
the relative strength of the CO rotational line profiles 
(jPecin et a l. 2006). 

Our modelling results suggest that the mass loss is 
varying as ^ 4.5 10~^ (10^^ cm/r)"-^ Mq/ji in the slowly 
expanding envelope (Vcxp = 15 km s~^), i.e. that the 
mass loss from VY CMa is increasing with time, reaching 
~ 6.4 10~^ Mq/jt in the inner part of the envelope (Rin) 
or about 100 years ago. In the bipolar lobes, the density 
enhancement between 5-7 10^^ cm represents an episode 
of strong mass loss about 350 to 500 years ago, for the 
expansion velocity of 45 km s~^. The nominal mass loss 
rate at 6 10^^ cm is 10"'^ Mq/jt. A similar episode of 
intense mass loss has also b een inferred from modelling 
results of iDecin et al.l (|2006[ ). although their assumption 
of spherical symmetry is inadequate to describe the ge- 
ometry and kinematics of the nebula around VY CMa. 

The current mass loss from VY CMa is better traced by 
dust continuum emission in the mid infrared or millime- 
ter and sub-millimeter wavelengths. We derived a mass 
loss rate of 4.7 10~^ Mo/yr from 1.3 mm dust continuum 
emission. Th i s valu e is very close to that obtained by 
iHarwit et al.l (|2001l ) from fitting infrared sp ectroscopic 
data. As already noted bv lDecin et al.l(|2006h . the differ- 
ence between mass loss rate derived from CO line profiles 
and dust continuum suggests that the dust continuum 
emission originates from the inner dense and hot region 
(a few tens of stellar radii from the central star), which 
is not well traced by CO low rotational lines. 

5.3. Origin of the complex envelope around VY CMa 



The origin of the bipolar outflow and the north-south 
elongated structure in VY CMa is still unclear at present. 
However, we note the similarity in morphology (bipo- 
lar outflows and/or disk- like structures) between the en- 
velope of VY CMa and that aro und several AGB and 
post- AGB star s such as tt^ Gru (IChiu et al.ll2006D. V 
Hya jHirano et al. 2004 ), M 2-56 (jCastro-Carrizo et all 
I200H) and X Her (NakashimaN2005l). The current popular 
model proposed by Morris (1987) to explain the striking 
bipolar morphology and the presence of coUimated fast 
outflow in AGB and post-AGB envelopes involves a low 
mass binary companion. The companion can have strong 
effect on the structure of the envelope. Its gravitational 
pull can focus the AGB wind toward the orbital plane to 
create an equatorial density enhancement. In addition, 
the accretion of wind material into the companion leads 
to the formation of an accretion disk, which might drive 
the coUimated high velocity outflow. Although the over- 
all morphology of the envelope is similar to that around 
other AGB and post-AGB stars, the high velocity wind 
from VY CMa, as suggested by our data and modelling 
results, is only weakly coUimated with very wide opening 
angle. It is conceivable that the interaction between an 
unseen companion and VY CMa could produce equato- 
rial density enhancement, possibly identified as the elon- 
gated structure and induce the coUimated fast wind to 
precess. Subsequent interaction with the nebula would 
result in wide bipolar lobe. The prominent curved neb- 
ulous tail seen in scattered light and the extended and 
curved structure seen in our CO (2-1) data might actu- 
ally trace the prec essing outfiow. To ex plain the promi- 
nent nebulous tail. lMonnier et ah! (|1999D also suggest the 
existence of a rotating jet and estimate the rotation pe- 
riod of a few thousands years. That time scale might 
be related to the orbital period of an unseen companion 
aro und VY C Ma. However, high resolution images from 
ISmith et al.l |2001[ ) suggest that this structure is actually 
much more complicat ed than expected from the model of 
iMonnier et al.l ()1999[ ) and probably more consistent with 
a localized ejection event on the stellar surface. Our cur- 
rrcnt SMA data cannot distinguish between these possi- 
bilities. Future high angular resolution observations and 
monitoring are necessary to determine the distribution 
of material and kinematics in order to test the presence 
of a companion. 

VY CMa might have magnetic activity, resembling 
solar fiares or eruptions, as suggested by ' Smith et al.l 
()2001h in order to explain the presence of numerous arcs 
and bright knots seen in scattered light in the nebula. 
We speculate further that the bipolar flow with wide 
opening angle and the equatorial density enhancement 
could be similar to what is observed in the solar wind. 
The fast solar wind originates from coronal holes located 
near the poles where magnetic field lines are open and 
the slower but denser wind comes from the equatorial re- 
gion. The geometry of the solar wind is closely coupled 
to the dipolar field of the Sun. In the case of VY CMa, 
the presence and strength of the magnetic field can only 
be inferred indirectly from observations of SiO, H2O and 
OH masers. Observations of Zeeman effect in SiO maser 
fines suggest a strength of - 2 G (|Herpin et al.ll200^ 
for the magnetic field aroun d VY CMa. U sing simple 
estimates, Soker. (1998) and ISoker fc ZoabU (|2Cl02) sug- 
gest that a magnetic field of about 1 G can significantly 
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influence the outflow. The observed strength of the mag- 
netic field in VY CMa thus suggests that it might affect 
the dynamics of the wind and if the field is organised on 
large scale, i.e. dipolar configuration, it could explain 
the structure of the envelope as seen in our SMA data. 

6. CONCLUSION 

In this paper, we present high angular resolution ob- 
servations of the circumstellar molecular envelope of the 
red supergiant VY CMa obtained with the SubMillime- 
tcr Array. Thanks to the wide frequency coverage of the 
SMA, we were able to trace the envelope in the ^^C0(2- 
1), ^'^C0(2-l) and 80(65-54) lines and dust continuum 
emission. We detected a strong and unresolved contin- 
uum emission at the center of the nebula. The derived 
dust mass is 1.5 10""^ Mq, implying a lower limit of the 
recent mass loss of 4.7 10"'* Mq/jt over the last 300 
years. The C0(2-l) emission is resolved into two distinct 
kinematic components: a slowly (15 km s~^) expanding 
shell elongated in the north-south direction and a high 
velocity (45 km s~^) component with a clear velocity gra- 
dient along the east-west direction. Similar structures 
are also seen in the ^^C0(2-l) and 80(65-54) emission. 
We interpret the high velocity component as a bipolar 
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outflow with a wide opening angle 120°), for which 
the symmetry axis is oriented close to the line of sight 
(i = 15°). Using a simple model, we could explain the 
main features of the 00(2-1) emission in our SMA data 
as well as the previous single-dish CO multi-line observa- 
tions. Density enhancement in the outflows, representing 
a previous episode of high mass loss, is inferred from our 
model. Our data complement previous high resolution 
optical spectroscopy and optical/lR images and provide 
a better understanding of the complicated morphology 
of the nebula around VY CMa. Future high resolution 
and high sensitivity observations of CO lines, including 
high-J transitions, will allow us to study the structure of 
the envelope and its physical properties in more detail. 
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Table 1: Basic data of VY CMa 

R.A. (J2000) 07''22™58!3315 (T)~ 

Dec. {J2000) -25°46'03yi74 (1) 

Distance 1.5 kpc (2) 

Angular scale 1" ~ 1500 AU ~ 2.25 10^^ cm 

Luminosity L, ~ 2 10^ Lq (3) 

Effective temperature Te//,, 2800 K (3) 

Stellar radius R* 1.2 lO" cm (3) 

Stellar mass M, > 15 Mq 

Mass loss rate M ~ 2 - 4 10—^ Mp/yr (4) 



Reference s. — (1) from the Hipparcos c atalogue , jPerryman ct al. (199?); (2) ILada fc ReidI II1978I '): (3) IMonnier et al.l l|1999l '): (4) 
[Bowers et all l[l983) ■ iDahchi et al.l l[T994 '). iStanek et al.l lfT995. ). ,Harwit et al., (,20011 ). 



Table 2: SMA line observations towards VY CMa. 



Line 


Rest freq. 
(GHz) 


Synthesized beam 
(arcsec) 


RA. 

(deg.) 


Vel. res. 

(km s-l) 


Icr rms 
(mjy/beam) 


1^00(2-1) 


230.538 


1.94 X 1.67 





4.2 


140 


i3CO(2-l) 


220.399 


2.65 X 2.30 


147 


11 


65 


SO(65-54) 


219.949 


1.99 X 1.80 


155 


4.4 


110 



Table 3: Parameters of our model 



Envelope inner radius (Rin) 5 10 cm 

Envelope outer radius (Rout) 8 10^^ cm 

Bipolar lobe inner radius (Riobc) 1-5 10^^ cm 

Opening angle of the bipolar lobes 120° 
Envelope expansion velocity 15 km s~^ 

Expansion velocity in the bipolar lobes [30 - 10 tan6/tan(60°)] (r - RiobG)/(Rout - Riobe) + 15 km s~^ ^ 
Gas density in the envelope 4.5 10'^ (10^*^ cm/r)^'^ cm~^ 

Gas temperature in the envelope 500 (2.25 10^^ cm/r)"'^ K 

Gas density in the bipolar lobes 6 10* cm~^ for Riobc < r < 5 10^® cm 

10^ cm"^ for 5 lO^*' cm < r < 7 10^® cm 

10^ (7 10^*^ cm/if (cm-^) for r > 7 10^® cm 

Gas temperature in the bipolar lobes 150 (1.5 10^® cm/r)"'^ K for Riobc < r < 5 10^'' cm 

35 K for 5 10^® cm < r < 7 10^^ cm 

35 (7 10^^ cm/r) K for 

Inclination angle 15° 

Position angle of the nebula axis 90° 

f=[C0]/[H2] 3 lO"-* 

^'^C/^'-^C 60 

Total molecular gas mass in the envelope 0.15 Mq 
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Fig. 1. — 1.3 mm continuum emission map of VY CMa. The synthesized beam is 1.38" x 1.23" (PA = 169°) as in dicated by the ellipse 
in the bottom left corner. Contour levels are every 40 mjy/beam {5a). The positions of the SiO HWright et al.llT990l) , H2O ||B owers et al.l 
[1993) and OH l IBowers et al.|[T983) maser centroids are indicated. The 1.3 mm continuum peak is located at R.A. = 07''22™58!339 and 
Dec. = -25°46'03"24 (J2000). 
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Fig. 2. — Channel maps of the ^■^00(2-1) line emission. The synthesized beam sizes 1.94" x 1.67" (PA = 0°) and is shown in the bottom 
left corner of the figure. Contour levels are drawn every 0.55 Jy/beam or 3.9 K (~ 4ct). The velocity resolution is 4.2 km s"'^. The cross 
indicates the position of the phase center. 
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Fig. 3. — Channel maps of the ^'^00(2-1) hnc emission. The synthesized beam sizes 2.65" x 2.30" (I'A = 117°) and is shown in the 
bottom left corner of the figure. Contour levels are drawn every 0.2 Jy/beam or 0.82 K (~ 3it). The velocity resolution is 11 km s~^. The 
cross indicates the position of the phase center. 
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Fig. 4. — Channel maps of the 80(65-54) line emission. The synthesized beam sizes 1.99" x 1.80" (PA = 155°) and is shown in the 
bottom left corner of the figure. Contour levels are drawn every 0.44 Jy/beam or 3.1 K (~ 4<t). The velocity resolution is 4.4 km s~^. The 
cross indicates the position of the phase center. 
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Fig. 5. — Integrated maps of the i2co{2-l), 13C0(2-1), 80(65-54) emission lines. The synthesized beam sizes are indicated in the 
bottom left corner of each frame. Contour levels are every 140 {7a), 30 (So-), and 105 {7a) K km s~i respectively. Offsets refer to the 
position of the phase center. 
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Fig. 6. — Position - velocity diagrams along the east-west and north-south directions for the ^■^CO(2-l) and 80(65-54) line emissions. 
Countour levels every 0.60 and 0.48 Jy/beam (4(t), respectively. 
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LSR velocity (km/s) 

Fig. 7.— Spectra of the ^^CO(2-l), ^=^0 0(2-1) and 80(65-54) lines smoothed to a 19'.'7 FWHM beam. The i2cO(2-l) spectrum observed 
with the JCMT bv lKemper et al.l (2003) is shown as a dashed line in the upper panel. The velocity resolution is 4.2, 11 and 4.4 km s ^ 
from top to bottom. 
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Fig. 8. — Sketch of the geometry of our model. The central star VY CMa (S) is embedded in a spherically synimetrie slowly expanding 
envelope (light grey). Part of the nebula is occupied by a bipolar outflow with a wide opening angle (dotted pattern). Regions of density 
enhancement within the bipolar lobes arc also included (filled bands). The line of sight OS intercepts the bipolar outflow axis with an 
inclination angle of 15°. A stellar photon is scattered at point P inside the envelope to the observer. The angle w is defined as OSP. 
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Fig. 9.— Spectra of the CO(2-l), CO(3-2), CO(6-5) and CO(7-6) lines observed bv lKemper et al.l pOOl l using JCMT are shown (full 
lines) together with the predicted CO line profiles from our model (dotted lines). 
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Fig. 10. — SMA spectrum of the ^^CO(2-l) line convolved with a 19'.'7 FWHM beam (full line) together with the output spectrum from 
our model (dotted line). The model spectrum is smoothed to the same velocity resolution of the SMA data (11 km s^^). A ratio ^^C/^^C 
= 60 was inferred from our modelling. 
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Predicted channel maps of the CO(2-l) emission from the envelope around VY CMa. Contour levels are the same as in Fig|2] 
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Fig. 12.— 
direction. 



E ^ W N ^ S 

Predicted position-velocity diagram of the CO(2-l) emission along the east-west (left frame) and north-south (right frame) 



